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Primordial black holes (PBHs) are theoretical black holes which can be formed during the ra¬ 
diation dominant era through the gravitational collapse of radiational overdensities. It has been 
well known that in the context of the structure formation in our Universe such collapsed objects, 
e.g., halos/galaxies, could be considered as bias tracers of underlying matter fluctuations and the 
halo/galaxy bias has been studied well. Applying such a biased tracer picture to PBH, we inves¬ 
tigate the large scale clustering of PBH and obtain an almost mass-independent constraint to the 
scenario that the dark matter (DM) consists of PBH. We focus on the case where the statistics of the 
primordial curvature perturbations is almost Gaussian, but with small local-type non-Gaussianity. 
If PBHs account for the DM abundance, such a large scale clustering of PBHs behaves as nothing 
but the matter isocurvature perturbation which is strictly constrained by the observations of cosmic 
microwave backgrounds (CMBs). From this constraint, we show that, in the case where a certain 
single field causes both CMB temperature perturbations and PBH formation, the PBH-DM scenario 
is excluded even with quite small local-type non-Gaussianity, |/nl| ~ 0(0.01). 


I. INTRODUCTION 

The origin of the dark matter (DM) is one of the most 
earnest questions for scientists. A lot of candidates have 
been proposed, but it still remains an open question. A 
primordial black hole (PBH) [1] is one of the DM can¬ 
didates, which has been studied for a long time. PBHs 
are black holes which are theoretically suggested to be 
formed in the early universe, and distinguished from or¬ 
dinary black holes which are believed to be formed from 
the collapse of stars. 

The constraint on the abundance of PBHs has been 
basically obtained from the non-detection of them (see 
Ref. [2]), and recently, by the Kepler telescope, all of the 
mass windows for PBHs to be a dominant component 
of DMs would be closed and the possibility of PBH-DM 
scenario seems to be completely excluded if PBHs have 
a monochromatic mass function [3, 4]. However, there 
remain theoretical uncertainty of constraints, the possi¬ 
bility such that the mass region of PBHs broadens widely, 
and so on. 

In this paper, we consider the large scale clustering 
behavior of PBHs and discuss the possibility of PBH- 
DM scenario in a way independent of the PBH mass. 
Here, we assume the case that PBHs are formed from 
the gravitational collapse of the radiational overdense re¬ 
gions, though there are several mechanisms for the PBH 
formation like a collapse of cosmic strings [5] and domain 
walls [6], or bubble collisions [7]. 

In the context of the structure formation, the large 
scale clustering of such collapsed objects as halos/galax¬ 
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ies, has been studied well. For example, Bardeen et al. [8] 
showed that the collapsed objects tend to be formed in¬ 
tensively in widely dense regions, in the peak-background 
split picture in which source density perturbations are 
divided into short- (peak) and long-wavelength (back¬ 
ground) modes. In the case where the statistics of the 
initial primordial perturbations is Gaussian, the large 
scale clustering of the collapsed objects could be related 
to the distribution of the underlying matter density field 
by introducing a constant scale-independent bias param¬ 
eter. On the other hand, due to the existence of the non- 
Gaussian feature in the statistics of the initial primordial 
perturbations, the bias parameter tends to have a scale- 
dependent component, which is known as scale-dependent 
bias, and this has been considered to be a powerful tool 
to hunt the primordial non-Gaussianity and extensively 
studied (see e.g., Ref. [9]). 

With the object that PBHs are also collapsed objects, 
Chisholm [10] considered the scale-independent bias of 
PBHs and concluded that though the bias factor is large 
in itself, the density perturbations of PBHs are hardly 
produced on much larger scales like observable scales 
of cosmic microwave backgrounds (CMB) temperature 
anisotropies. Recent several papers [11, 12] also agree 
with this result. This is simply because that PBHs are 
caused by the collapse of horizon scale overdensity, which 
are hardly affected by the super-horizon density fluctu¬ 
ations that has little causality with physics inside the 
horizon. 

In this paper, based on these previous works, we con¬ 
sider the case where the statistics of the primordial curva¬ 
ture perturbation is almost Gaussian, but with small non- 
Gaussian component, and investigate the large scale clus¬ 
tering of the PBHs. As a non-Gaussian component, we 
consider the so-called local-type non-Gaussianity, which 
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is known to give a strong scale-dependence in the halo 
bias. We find that if PBH is a dominant component of 
DM, the large scale clustering of the PBH can be de¬ 
tected as the matter isocurvature perturbations and ex¬ 
cluded by CMB observations even in the case with small 
non-Gaussianity. 

This paper is organized as follows. In section II, we 
briefly review the bias effect and calculate the matter 
isocurvature perturbations from biasing effect of PBH in 
the case where a single source causes both observed CMB 
temperature perturbations and PBH formation. In sec¬ 
tion III, we discuss the other case where the source of 
PBH formation is different from that of CMB perturba¬ 
tions. Section IV is devoted to the conclusions. 


II. SINGLE SOURCE CASE 

In this section we discuss the simplest case that the 
fluctuation of a certain single field causes both observed 
large scale CMB temperature perturbations and PBH 
formation. Let us denote a primordial adiabatic cur¬ 
vature perturbation on the comoving slice as TZ. Al¬ 
though the amplitude of 1Z is known to be around 1CV 5 
on CMB scale, it could be large enough to produce PBHs 
on smaller scales, becoming blue-tilted [13, 14] or having 
some bump on its power spectrum [15-18]. Also we as¬ 
sume TZ has a small local-type non-Gaussianity which is 
represented as 

^(x) = 5 (x) + /nl(s 2 (x) - (g 2 )), (1) 

where g is a Gaussian field. Here, we do not mention 
the concrete PBH formation model where the generated 
adiabatic curvature perturbation has the above type of 
non-Gaussianity. There have been several claims about 
the local-type non-Gaussianity in single source case. Ac¬ 
cording to Maldacena [19], there is a consistency relation 
between non-linearity parameter /nl and the spectral 
tilt of the adiabatic curvature perturbations, and hence 
even in the simplest inflation model, namely single-field 
slow-roll inflation, the local-type non-Gaussianity does 
not vanish and is O(10 -2 ). On the other hand, recently, 
some authors claimed that such a non-Gaussianity given 
by the consistency relation is nothing but a gauge artifact 
and it does not contribute the scale-dependent bias [20- 
22] and this statement is based on the adiabaticity of 
the curvature perturbations. As for the PBH formation 
model, e.g., in the bump cases, in order to grow up the 
curvature perturbations rapidly on the bump scale the 
adiabaticity would be violated and it would be possible 
to generate the curvature perturbations with non-zero 
local-type non-Gaussianity. Related issues in preheating 
model was discussed in Ref. [23] . Anyway we simply pre¬ 
sume a small non-Gaussianity /nl ~ 0(0.01) here with¬ 
out any assumption of the concrete inflationary models. 


A. PBH bias 

We will assume that PBH formation can be described 
by Press-Schechter approach [24] as in the standard 
structure formation. Following the peak background split 
picture in the structure formation [8] , we discuss the large 
scale clustering of PBH. In other words, we will con¬ 
sider how the PBH number density in some region of 
scale R will be modulated if that region is on some long- 
wavelength fluctuations whose scale Ri is much larger 
than R. Anyway let us calculate the PBH density with¬ 
out the long-wavelength mode at first. We would like to 
note that using the peak background split picture is just 
for intuitive understanding how the non-Gaussianity af¬ 
fects the large scale clustering of the biased objects such 
as halos, galaxies and PBHs. A clustering behavior of 
halos/galaxies with the non-Gaussian initial fluctuations 
has been also discussed by employing not only the peak- 
background split picture but also other analytic approach 
[25, 26]. There are also several literatures about numeri¬ 
cal N-body simulation with the non-Gaussian initial fluc¬ 
tuations and they are consistent with the analytic result 
(see, e.g., Ref. [9]). 

According to [11], when the PBH formation is dis¬ 
cussed, one should not use the curvature perturbations 1Z 
but take the density perturbations of the radiation on the 
comoving slice 5. The Fourier transformed component of 
such density perturbations of the radiation is obtained 
from the Poisson equation and it becomes [11] 

i < k > = K^) 2R<k) - (2) 

on super-horizon scale. Following the Press-Schechter 
approach, we consider the density perturbations coarse¬ 
grained on PBH-formation scale. With some window 
function W(kR s ), we define the coarse-grained density 
perturbations as 

d s (k) =M s (k)lZ(k), (3) 

where 

M a {k) = \ (—p- ) W(kR s ) = \(kR PBH ) 2 W(kR s ). (4) 

9 \aH Ipbh/ 9 

R s is the coarse-graining scale, and almost the same as 
the horizon scale at the PBH-formation time Rpbh = 
(aiijpBH) -1 since the overdensity will collapse soon after 
entering the horizon. Note that R should be much larger 
than R s . 

In the Press-Schechter approach, the overdensity such 
that <5 s (x) exceeds a certain threshold 8 C is assumed to 
collapse to be PBHs. 1 Therefore, assuming the statistics 


1 In the case such that 1Z has the non-Gaussianity, the dynamics 
of collapse after the overdensity enters the horizon may also vary. 
Therefore we may not be able to say sweepingly the overdensity 
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of S s is also nearly Gaussian, 2 the probability of PBH 
formation is given by 


Pi(> v) = 


y/2r 


diexpl-^ 


= 6rfC ,V2 ~ 




( 5 ) 


where v = 5 c /a s and cr s is the standard deviation of S s 
averaged in the region R: 

v 2 s = ( S s)r = y J R d 3 xS s ( x ) 

= y J R d 3 x I ^l Ms {k)M s {q)H^)n{ qK^)- 

( 6 ) 


where V is the volume of the region R. We have in¬ 
cluded conventional factor 2 of Press-Schechter formula¬ 
tion. Also we have assumed a high peak limit as v 1, 
because PBHs should be rare objects. Then we will con¬ 
sider how this probability will be modulated under the 
existence of long-wavelength modes. 

With the long-wavelength mode, the two quantities S c 
and ct s seem to be able to be modified. However, in the 
integral of Eq. (6), the factor ( k/ciH ) 2 in M s {k) cuts off 
the long-wavelength modes of 7Z and er s is determined 
only by the short-wavelength modes as 


a 2 = 


r- A R-l 


d\ogkVn(k)M 2 s (k)~V n (k s )M 2 s (k s ),(7) 


with the power spectrum Viz(k ) = |U(k)ft(-k) = 
^\lZ(k)\ 2 . Here we have replaced V~ x f R d 3 a;e i ( k+q > x 
by (27r) 3 i5^ 3 ^ (k + q) for sufficiently large k and q com¬ 
pared to R~ 1 . If 1Z is completely Gaussian and there is 
no modal coupling between long- and short-wavelength 
modes, it does not affected by long modes at all. On the 
other hand, the threshold S c is simply reduced to S c — Si 
if the region R is on the long-wavelength density pertur¬ 
bation Si (see Fig. 1). Therefore the PBH density in that 
region will be slightly larger as 


< 5 pbh(x) 


Pi(> v\Si(x.)) 


Pi(>v) 

dS c d log Pi 

dSi d5 c 


M x ) 


1 


V 

Os 


Si(x)- 


( 8 ) 


such that 5 a > 5 C will collapse (For instance, Nakama et al. [27] 
analyzed the development of overdensities of generalized curva¬ 
ture profile beyond the Gaussian one and proposed two crucial 
parameters as thresholds of PBH formation.). However those ef¬ 
fects are beyond the scope of this paper, and we just follow the 
standard Press-Schechter approach. 

2 The effect of the non-Gaussian profile of 8 a is considered in [28— 
31] for example. However such an effect to the scale-dependent 
bias is the second order quantity w.r.t. f\ i [32-34], and there¬ 
fore it can be neglected safely. Moreover, in appendix A, it is 
suggested that the part of bias factor without f\ i is not so dif¬ 
ferent even if the non-Gaussian profile effect is taken into account 
as long as the PBH abundance is fixed. 
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FIG. 1. Schematic diagrams of a bias effect in the peak- 
background split formulation. Gray dotted lines represent 
the threshold for PBH formation S c and dark blue points show 
the regions which will be PBHs. In the first (top) figure, it 
can be seen that the threshold for the short-wavelength mode 
S a is effectively reduced to 8 C — Si. Therefore PBHs tend 
to be formed in Si > 0 region, and this effect is called the 
scale-independent bias. However, long-wavelength modes of 
the comoving radiational perturbations are practically quite 
suppressed compared to the curvature perturbation IZi due 
to the factor (kRpBn) 2 in Eq. (4) as the second (bottom) 
figure shows. Even in such a case, if /nl > 0 (or < 0), the 
amplitude of S a itself tends to be larger, being proportional 
direct to 7 Zi via the non-Gaussianity of primordial curvature 
perturbations. Such a bias effect is called the scale-dependent 
bias. 


Here we have used a high peak limit v 1 again. The 
spatial label x denotes the coarse-grained position of the 
region R. This coefficient b 0 := — d Pc ~ is called 
scale-independent bias. 

Thus, following the peak-background split picture, we 
find that the number density fluctuations of PBH, < 5 pbh, 
on large scales is given by <5 pbh — jf-S in the pure Gaus¬ 
sian case. However, although a biasing factor v/a s is 
larger than order of unity, dpBH is actually much smaller 
than 7 Z on the CMB scale. This is because the CMB scale 
^CMB Quite larger than the PBH scale -Rpbhj namely 
fccMB-RpBH 1, and then (5 pbh is strongly suppressed 

as <5 pbh(^CMb) ~ v/a s X (&CMBl?PBH) 2 7'l(fccMB) ^ 

TZ(kcMB)- Indeed, this result is consistent with the 
claims of [10-12]. However, it will be dramatically 
changed if we consider the effect of non-Gaussianity, 
namely the scale-dependent bias [9]. 

If the primordial perturbations have local-type non- 
Gaussianity as given in Eq. (1), there is a correlation be¬ 
tween short- and long-wavelength perturbations as shown 
later. Therefore the long-wavelength density perturba¬ 
tion Si not only reduces the threshold effectively but also 
modifies the variance of short-wavelength perturbations, 
a s . Thus the bias parameter b = d 1 q| Pl ^ has the 
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other component 

( 9 ) 

Si=0 

which is called scale-dependent bias. Then let us evaluate 
this bias effect. 

From Eq. (1), TZ(x) can be decomposed into long- and 
short-wavelength components as 

H(x)=gi (x) +g s (x) +/nl [(gi (x) +g s (x)) 2 - (g?) - (g 2 )} 
= Sz(x)+/nl(s?(x)-(& 2 )) 

+9s (x) + /nl (2 gi (x)g s (x) +g 2 s (x) - (g 2 )). (10) 


_ d log <r s d log Pi 
dSi d logcr s 


In the second equation, the first line which is independent 
of g s denotes the long-wavelength mode TZi and the sec¬ 
ond line is regarded as the short-wavelength mode TZ S . 
Til will be cut off in the integral for a 2 (6) due to the 
transfer function A4 s (fc) as we have mentioned, but now 
we have another contribution of long-wavelength modes 
to the variance via gi(x) in 1Z S . Since the second or¬ 
der terms of g s (x) do not include the long-wavelength 
mode g/(x) and their contributions to bias effects are 
suppressed as the second order of /nl [32-34], one can 
neglect these terms and see that the variance a 2 in the 
region R is just constantly amplified as 

ct s (x) = (1 + 2 / NLi 9 z(x))CT s , (11) 


noting that gi(x) is almost constant in the region R. Here 
d 2 denotes the variance of S s averaged over the whole uni¬ 
verse, namely d 2 = ( S 2 ). Hence, with an approximation 
Til (x) ~ 9 /(x), the Fourier mode of the dependence of a s 
to 5i on large scale is given by 


Slog a s (k) 

dSfik) 


6 z ( fe )=0 


Mi\k) 


Slog a s {k) 
dgi{k) 


gi(k)=0 


2h L MY\k), (12) 


where, 


M t (k) = ^(kR. PBU ) 2 W(kRi). (13) 


On the other hand, since Pi is the function only of v, we 
can obtain the following relation with use of the definition 
of b 0 . 


Slog Pi dv dlog Pi 
S log <t s s dcr s dv 


S c 


S log Pi 

~dS7~ 


S c b 0 . (14) 


From the above results, the scale-dependent bias is given 
by 

Ab(k) = 2f NL MT 1 (k)6 c b 0 . (15) 


It is noteworthy that A b has a factor of yVI / 1 . Since 
this factor cancels out (fcPpBH) 2 in S s , (5 pbh with the 
scale-dependent bias is not negligible compared to the 
adiabatic curvature perturbation 1Z even on much larger 
scales. This is because the scale-dependent bias is di¬ 
rectly proportional to TZ via the non-Gaussianity, while 
the scale-independent bias is proportional to S. 


B. PBH-DM isocurvature perturbation 

Gravitational collapse is irreversible process and po¬ 
tentially accompanied by entropy production. Further¬ 
more, especially in the case of PBH formation, the equa¬ 
tion of state of the density component changes from the 
radiation to the matter, and therefore the density per¬ 
turbations of PBHs can bring about easily-detectable en- 
tropic (or isocurvature) perturbations. In this section we 
will calculate the specific value of the isocurvature per¬ 
turbations. 

If a dominant component of DMs consists of PBH, the 
matter isocurvature perturbation can be given in terms 
of 5 pbh as 

S := <5pbh — ^S. (16) 

By using Eq. (15), the Fourier transformed component 
of the matter isocurvature perturbation smoothed on the 
scale R is 

S s (k)=(t 0 +Ai(t) —j)fe(k) 

= (b 0 M R (k) + 2f NL S c b 0 ^(k)-^M R (k^7Z(k), 

(17) 

where M R {k) = ^(kR P Bn) 2 W(kR). As mentioned be¬ 
fore, the factor (fccMB-RpBH) 2 br A4 _r is so small that we 
can neglect the first and third terms in the parenthesis. 
Thus the power spectrum of the isocurvature mode is 
given by 3 

Ps(kcMB) — (2f^hS c bo) 2 Pn(kcMB) 

— (2/nl^ 2 ) Pr(^cmb)- (18) 

Here we have approximated ^y(fccMB) as unity be¬ 
cause is larger than both Ri and R now and then 

W(kcMBRi) — W(kcMBR) — 1- Note that S R is sim¬ 
ilarly equal to the unsmoothed mode S on CMB scale. 
Also we have used v = 5 c /a s and bo ~ v/a s in the second 
equation. 

CMB observations have already given tight constraints 
on the amplitude of the matter isocurvature perturba¬ 
tions and, according to the Planck collaboration [35], in 
the case of fully (anti-) correlated type which corresponds 
to S oc TZ (or S oc —TZ), the matter isocurvature mode is 
constrained as 

Vs < f 0-0025 (fully correlated, /nl > 0) 
Vs+Vn ~ | 0.0087 (anti-correlated, /nl < 0), 


3 Note that this result does not mean the rarer PBHs bring about 
the more isocurvature perturbations even though the higher v 
means the lower abundance of PBHs. This result supposes 
flpBH — ^DMi an d if PBHs are sub-component of DMs, the 

isocurvature modes are proportional to ( ) w ^ 1 ^ c ^ 1 

fact exponentially suppressed as shown in Eq. (5). 
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at 95% CL. Therefore, from Eq. (18) we must satisfy 


|/nl|^ 2 ~ 


VO.0025 

2 

V0.0087 

2 


0.025. (/nl > 0) 
0.047. (/nl < 0) 


( 20 ) 


Assuming the nearly monochromatic mass function for 
PBHs, the current abundance of PBHs is given by [2] 


/ m \ -1 ^ 2 

L!pbh ~ 0.86 x 10 8 Pi i—j , (21) 

where Mq denotes the solar mass ~ 2 x 10 33 g. Inversely, 
if DMs consist of PBHs and Hpbh = Hdm = 0.31 [36], 
the probability Pi should satisfy the following relation 

/ m \ 1 ^ 2 

P! ~ 0.36 x 1CT 8 (—) ■ (22) 


Furthermore, from Eq. (5), v can be written in terms of 
Pi as 

v = V2erfc _1 (Pi), (23) 

where erfc -1 denotes the inverse function of the com¬ 
plementary error function. Therefore, if we assume the 
mass of PBH is less than 10 4 ° g, 4 we obtain the following 
relation 


V>2 


erfc -1 0.36 x 10 -8 


10 4 ° \ 
2 x 10 33 / 


1/2N 


~20. (24) 


Here, note that erfc (x) is a monotonic decreasing func¬ 
tion. It obviously conflicts with the constraint (20) with 
/nl ^ 0(0.01). 

As a result, we find that it is difficult to satisfy 
the isocurvature constraint (20) even with a small non¬ 
linearity parameter |/nl| ~ 0(0.01) and this means that 
the PBH-DM scenario is strictly constrained due to the 
existence of the small local-type non-Gaussianity in the 
single sourced case. This is a main result of this paper. 

Furthermore, we have also considered the case that 
PBH is a subdominant component of DMs. We have 
plotted the upper limit of the PBH fraction in such a 
case to satisfy the isocurvature constraints (19) in Fig. 2. 
Note that the isocurvature power and the probability of 

PBH formation are given by Vs = ( Qd’m ) ^ pbh an d 
Pi ~ 0.36 x 10- 8 ^ J respectively in that case. 


4 Here we respect the plot of Ref. [4] in which the PBH constraints 
were summarized up to Mpbh = 10 40 g. Moreover Afshordi 
et al. showed that for more massive PBHs than ~ 10 37 g, the 
Poisson-noise fluctuations deviating from the adiabatic pertur¬ 
bations become too large to be consistent with observed ag [37]. 
Therefore we constrain the PBH mass to be less than around 
these values. Also, lighter PBHs than 10 15 g have been already 
evaporated due to the Hawking radiation by today, and therefore 
we are interested only in more massive PBHs than 10 15 g 
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FIG. 2. The plot of the upper limit of the current PBH 
fraction Dpbh/Ddm from the isocurvature constraints (19). 
Shaded regions represent existing constraints from various ob¬ 
servations [3, 4]. The black solid and dotted lines show the 
upper limit in the case of /nl = ±0.01 and the red solid and 
dotted ones exhibit the constraints in the case of /nl = ±0.1. 
Also from this plot, it can be seen that PBHs cannot be the 
main component of DMs even with small non-linearity pa¬ 
rameter |/nl| = 0.01. 


III. DIFFERENT SOURCE CASE 

In the case where the source fluctuation of PBH for¬ 
mation is different from that of CMB fluctuations, the 
above discussion does not change so much. Let 7Z and 
E denote the curvature perturbations which respectively 
induce the CMB fluctuations and the PBH formation. In 
other words, 1Z is dominant and has a nearly flat spec¬ 
trum on CMB scale, while E dominates on much smaller 
scale. 

If we assume E is almost Gaussian but has a small 
local-type non-Gaussianity as 

E(x) = g(x) + /nl,e(5 2 (x) - ( g 2 )), (25) 

the power of the matter isocurvature perturbations is 
given by 

Vs- (2/ N L,E p2 ) 2 iPs, (26) 

like as Eq. (18). By introducing an effective non-linearity 
parameter defined as 



the power of the isocurvature perturbation can be written 
as 

v s = (2/ N L,eff p2 ) 2 P^. (28) 

From the above expression, we can also obtain the con¬ 
straints on the PBH-DM scenario in the different source 
case depending on the ratio between the power of 1Z and 
that of E on CMB scale. Note that in this case the cor¬ 
relation between adiabatic curvature perturbations and 
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generated isocurvature perturbations is much suppressed, 
that is, almost uncorrelated-type isocurvature perturba¬ 
tions. The constraint for the uncorrelated isocurvature 
perturbations shown by Planck collaboration is [35] 

^’ S 0.036, (uncorrelated), (29) 

Vs + Viz 

at 95% CL. and therefore the constraint for /NL,eff is 
weakened up to 

|/NL,effk 2 ^ 0.095. (30) 

Furthermore, E is not necessarily almost Gaussian but 
possibly quite non-Gaussian. Accordingly, we also dis¬ 
cuss the chi-squared case in appendix A as a simple ex¬ 
ample. 


IV. CONCLUSIONS 


function. We leave more precise discussion about the 
effect of non-monochromaticity for a future issue. 

We also think that it should be interesting to evaluate 
the specific value of /nl in some concrete cases where 
abundant PBHs are produced. Further, though we fo¬ 
cused on simple local-type non-Gaussianity here, it will 
be interesting and important to discuss the bias effect 
for other types of non-Gaussianity. We also leave such 
studies for future issues. 
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In this paper, we consider the clustering of the PBHs 
on CMB scale produced by the bias effect. Since there 
is a great gap between CMB scale and the horizon scale 
at the time of PBH formation during the radiation domi¬ 
nant phase, the distribution of PBH is hardly biased basi¬ 
cally. However, if the curvature perturbation field which 
cause PBHs has a local-type non-Gaussianity, even if it is 
small, we have showed non-negligible spatial fluctuations 
of PBH number density can be produced on CMB scale 
because the source field has a correlation between large 
and small scale via the non-Gaussianity. Specifically, the 
fluctuation of the number density of PBH is biased from 
the source curvature perturbation roughly by the fac¬ 
tor of non-linearity parameter /nl times (erfc - (Pi)) , 
where Pi is the PBH production probability which is 
commonly represented as /3 in the literature. If DMs 
consist of PBHs, these extra density perturbations on 
CMB scale can be detected as the matter isocurvature 
perturbation and are constrained strictly by CMB obser¬ 
vations. 

We have showed, if the source curvature perturba¬ 
tion of PBHs also cause CMB temperature fluctua¬ 
tions, the PBH-DM scenario is excluded even with quite 
small non-linearity parameter |/nl| ~ 0(0.01). On 
the other hand, in the case that the source curva¬ 
ture perturbation of PBHs E is different from that of 
CMB fluctuations 1Z, we can avoid the isocurvature con¬ 
straints if E is sub-dominant enough on CMB scale. 
Even such a case, the PBH-DM scenario should satisfy 

l/NL,s|^^;(fecMB) ~ 0(0.001). These constraints are 
almost independent of the PBH mass, and moreover they 
may not change a lot even if the PBH-mass spectrum is 
quite non-monochromatic unlike many other constraints. 
This is because even if the isocurvature perturbations of 
each mass range are small enough, they will be piled up 
and the total amplitude of them will be determined only 
by /nl and Hpbh 7 hardly affected by the shape of mass 


Appendix A: Chi-squared type 


In this appendix, we consider the simple chi-squared E 
which is characterized as 


E(x) = ff 2 (x) — (g 2 ), (Al) 


as an example of fully non-Gaussian case. Suppose the 
coarse-grained density perturbations 5 S also follow the 
chi-squared distribution, the probability of PBH forma¬ 
tion is given by 


,oo 1 

P 1 (>5 C ) = 2 dS - 

“W yvr <j a (y/ 


( \/2S+a s ) 


: exp 


= 2 erfc 


' \[2v + 1 


^ 2<t s j 


— 2 * 


tt(\/2v + 1) 


exp 


y/2is + 1 


(A2) 


where v = 5 c /a s again. The scale-independent bias can 
be calculated similarly to section II as 


b Q = — 


d log P\ 

~dS7~ 


1 d log Pi 
a s 


1 

\/2 a s 


(A3) 


On the other hand, to calculate the scale-dependent part 
of the bias parameter, we have to discuss the dependence 
of (J s on Si. From Eq. (Al), we decompose E into long- 
and short-wavelength components like as Eq. (10), 

S(x)= ft 2 (x)-<s?> 

+^( x )-(ffs) + 2 5i( x )Ss( x )- ( A4 ) 


Similarly to the discussion in section II, the second 
line terms in the expression of E(x) denote the short- 
wavelength mode of E. Noting that third moments of 















7 


g s should vanish, the power spectrum of S s on i?-scale is 
given by 


^(fc)lE,(x)=^aW[ 4 (3?(x)-(ft 1 2 ))iP 5s (fc)+P Ss (fc)| S! (x)=o] 
= Ml(k) [4E, (x)P 9s (k) +7Y (k) | Sl(x) =o] ■ (AS) 

Then the standard deviation of S Sl which is the square 
root of the integration of Vs s , can be approximated by 


( 2er 2 

cr s (x) ~ cr s |s z ( x )=o I 1 + -^E/(x) 


(A6) 


where 


Therefore 


Vg,s = J dlogkV g (k)Ml(k). (A7) 


d log a s (k ) 


dStik) 


<5i=0 




(A8) 


and with use of Eq. (14), we obtain the scale-dependent 
bias as 


2cr 2 . 

A b= -77 


M~T 1 5 c b 0 ~ 


V^cr 2 


-M 


-1, 


(A9) 


Since it has the factor of M^ 1 again, the number density 
fluctuation of the PBH can be produced through this 
scale-dependent bias even on CMB scale and the power 
of matter isocurvature is given by 


'Ps(fccMB) 



2 


2 

'Ps(kcMB)- 


(A10) 


Introducing the effective non-linearity parameter as 

y/2cr 2 /- 

/NL,eff = fj 2 Jl^(fccMB), we obtain the constraint 


as 


/NL.effi' £ 0.095. (All) 


Note that v is not squared unlike the case of small non- 
Gaussianity and this fact just comes from the apparent 
difference of Pi, Eq. (5) and (A2). The lower bound for 
v given by 


v = 


2[erfc- 1 (P 1 /2)] 2 —1 

72 


(A12) 


is not so different from that of v 2 in the case of local-type 
non-Gaussian, Eq. (24). 

We can not calculate /NL,eff without determining the 
detail profile of E (or g ) like a spectrum index, but can 
understand this result qualitatively. First, let us approx¬ 
imate by ^(Rj 1 ). In fact, the contribution around 
smoothing scale R s is dominant for the integrations for 
a 2 s and a 2 due to the factor of M s - Moreover Vy can 
be approximated by V 2 . Therefore <J 2 s /cr 2 is roughly 
given by V ~ 1 . On the other hand, in the case of nearly 
scale-independent g , local non-linearity parameter /nl,e 
for E can be written as -Bs(fci,fc 2 ,fc 3 )/(PE(fci)Ps(fc 2 ) + 
2 perms.), where (27r) 3 4 5 6 7 8 9 <5 (3) (k 1 + k 2 + k 3 )P s (fc 1 , k 2 , k 3 ) = 
(E(ki)E(k 2 )E(k 3 )) and Ps(fc) = |E(k)| 2 . According 
to Wick’s theorem, we can obtain ~ 

P g (P g (fci)P 9 (fc 2 )+2 perms.). Moreover (Pg(ki)Pg(k 2 ) + 
2 perms.)/(Ps(fci)Ps(fc 2 ) + 2 perms.) ~ V ~ 2 if g is al¬ 
most scale-independent. Therefore cr 2 s /cr 2 ~ V ~ 1 in¬ 
deed denotes the effective local non-linearity parameter 
/nl,e in the case of chi-squared field. 
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